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> ABSTRACT 

CD We present a sample of 29 galaxy clusters from the XMM-LSS survey over an area of some 

5 deg^ out to a redshift of z = 1.05. The sample clusters, which represent about half of the 
-ray clusters identified in the region, follow well defined X-ray selection criteria and are all 
spectroscopically confirmed. For all clusters, we provide X-ray luminosities and temperatures 
as well as masses, obtained from dedicated spatial and spectral fitting. The cluster distribution 
peaks around z = 0.3 and 7=1.5 keV, half of the objects being groups with a temperature 
below 2 keV. Our Lx - T(z) relation points toward self-similar evolution, but does not exclude 
other physically plausible models. Assuming that cluster scaling laws follow self-similar evo- 
^ lution, our number density estimates up to z=l are compatible with the predictions of the con- 

• fH cordance cosmology and with the findings of previous ROSAT surveys. Our well monitored 

selection function allowed us to demonstrate that the inclusion of selection eff'ects is essential 
H for the correct determination of the evolution of the Lx - T relation, which may explain the 

contradictory results from previous studies. Extensive simulations show that extending the 
survey area to 10 deg^ has the potential to exclude the non-evolution hypothesis, but that con- 
straints on more refined ICM models will probably be limited by the large intrinsic dispersion 
of the Lx - T relation, whatever the sample size. We further demonstrate that increasing the 
dispersion in the scaling laws increases the number of detectable clusters, hence generating 
further degeneracy [in addition to erg, O^, - T(z)] in the cosmological interpretation of the 
cluster number counts. We provide useful empirical formulae for the cluster mass - flux and 
mass - count-rate relations as well as a comparison between the XMM-LSS mass sensitivity 
and that of forthcoming SZ surveys. 

Key words: surveys - X-rays: galaxies: clusters - large-scale structure of Universe - cosmo- 
logical parameters 



*Based on data collected with XMM, VLT, Magellan, NTT, and 
CFH telescopes; ESO programme numbers are: 070.A-0283, 070.A-907 
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1 INTRODUCTION 

Along with Cosmic Microwave Background (CMB) measurements 
and supernova (SN) observations, clusters of galaxies provide key 
cosmological information. It is especially instructive to cross-check 
the constraints from these three classes of data, since they originate 
from different physical processes. Moreover, since theory and nu- 
merical simulations allow us to follow cluster formation from the 
initial power spectrum, which is directly measured from the CMB, 
it is critical to test that the "CMB WMAP concordance cosmol- 
ogy" is consistent with the observed properties of clusters in the 
low-z Universe. 

In the framework of hierarchical cosmic structure formation 
involving Cold Dark Matter type scenarios (CDM), where the 
smallest perturbations collapse first, clusters correspond to the mass 
scale that entered the non-linear regime between redshift three and 
the present epoch. In this sense, the most massive galaxy clusters 
in the local Universe represent the largest virialised structures. This 
property of being both "relaxed" and rare 3cr events regime has 
been extensively exploited through formalisms like that of Press & 
Schechter (1974) in connection with the spherical collapse model, 
for the general case of Gaussian random field fluctuations. This 
connects, in an analytically tractable manner and for any redshift, 
both the cluster abundance as a function of mass and the cluster 
spatial distribution, to the properties of the initial fluctuation spec- 
trum - in particular its normalisation, erg, its shape, F, as a function 
of the density of the Universe, fl^. and the equation of state of Dark 
Energy (e.g., |Refregier et al|2002l|Majumdar & Mohr|2003l ). This 
first-order approach is well supported by numerical CDM simula- 
tions: clusters lie at nodes of the cosmic network, have virialized 
cores, and are still growing by accretion along filaments at a rate 
that depends on the cosmology. However, at the same time as clus- 
ters started being used as cosmological tools, it was realized that 
the interpretation of their observed abundance as a function of time 
is actually very much dependent on the evolution of the observable 
cluster properties themselves. In order to break this latent degener- 
acy in such a way that clusters can effectively be used as cosmo- 
logical candles, it is essential to understand how cluster properties 
impact on their delectability at any epoch. 

While galaxies constitute only a few per cent of the total clus- 
ter mass, about 80% of the baryonic mass resides in the tenuous 
X-ray emitting intra-cluster medium (ICM), settled in the cluster 
gravitational potential. Because cluster X-ray emission is extended, 
clusters are readily identified among the high-galactic latitude X- 
ray population, which is dominated by point-like AGN. However, 
whatever the detection method - optical or X-ray - the fundamen- 
tal question of how to relate observable quantities to cluster masses 
remains. This is crucial because it is the cluster masses that enter 
the theory of structure formation as generated by theoretical cos- 
mological calculations. This issue becomes particularly important 
outside the local Universe, since the evolution of cluster properties 
is not well established. 

The Einstein Medium Sensitivity Survey (EMSS), followed 
by REFLEX based on the ROSAT All-Sky Survey, and a number of 
serendipitous clusters surveys (RDCS, SHARC, MACS, \60deg^ 
etc) from ROSAT deep pointings, have provided the first 'cosmo- 
logical cluster samples' (see the synoptic plot in P ierre et al.| ( |2004| ) 
and references therein). The mass-observable relations used in the 
analysis of data from these surveys relied on the assumption of hy- 
drostatic equilibrium and on the (mostly local) observed Lx - T 
relation. The cluster selection function - a key ingredient - was 
modelled using a variable flux limit across the survey area. Under 



these hypotheses, the derived cosmological constraints appeared to 
be in agreement with the concordance model [see review in lRosati,"] 
^Borgani & Norman (2002 )] . 

Within the same period, deep GINGA, ROSAT and ASCA ob- 
servations of nearby clusters revealed that the Lx-T relation is sig- 
nificantly steeper than expected from purely gravitational heating 
( Arnaud & Evrard 1999 ), hence suggesting the presence of other 
heating/cooling sources such as feedback from star formation or 
AGN, in addition to the effects of cooling flows. This particularly 
affects the low-mass end of the cluster population - groups with 
temperatures [0.5 - 2] keV where the gravitational binding ener- 
gies are low. 

Nearby cluster observations at high spatial resolution by 
Chandra have also shown that the ICM is not the well-relaxed 
medium previously assumed: shocks, cold fronts, and bulk veloc- 
ities are seen even in apparently relaxed clusters (e.g. Mazzotta,| 
[Fusco-Femiano & Vikhlinin|2002|[Dupke & Bregman|26o6| ). High 
spectral resolution XMM pointings have led to a totally new ver- 
sion of the putative "cooling flow" scenario, where episodic heat- 
ing/accretion by the central AGN could play a key role in the ICM, 
preventing any central "cooling catastrophe". Measurements of the 
Lx-T relation for distant (z > 0.5) massive clusters are in progress, 
whether the evolution of scaling laws follows simple self-similar 
expectations is still hotly debated [see review in |Arnaud| ( [2005a| )]. 
All these results currently pertain to the upper end of the cluster 
mass function. They present a new challenge for high-resolution 
numerical simulations, which, in turn, should quantify deviations 
from hydrostatic equilibrium. 

In parallel, the building of large serendipitous XMM and 
Chandra cluster samples - from public archive data - has been initi- 
ated by a number of groups: the XMM cluster survey (XCS, ^Romer| 
|eran[200T]), SEXCLA S ( |Kolokotronis"erari|200gl ) and ChaMP 
( [Barkhouse et al. f 2006V Preliminary results over areas < 10 deg^, 
give cluster densities of the order 5deg~^ for objects detected in- 
dependently in both X-ray and optical wavebands. No cosmologi- 
cal analysis has been performed on these samples so far, but these 
searches have enabled the detection of the most distant X-ray clus- 
ters to date at z = 1.4 ( |Mullis et al.|2005| ) and z = 1.45 ( [Stanford et| 
|al.|2006| ). 

In addition to these large surveys, a number of contiguous 
surveys, such as COSMOS ( [Finoguenov et al.|2007| ) and XBootes 
( |Kenter et al.|2005l ), are also being conducted. Following on from 
the REFLEX cluster survey (Bohringer et al. 2002), and exploit- 
ing the unrivalled sensitivity of the XMM-Newton X-ray observa- 
tory, the XMM wide area survey [XMM-LSS, [Pierre et al.] ( [2004l )] 
is the largest contiguous X-ray cluster survey being undertaken at 
the present time. It has been designed to investigate the large scale 
structure of the Universe as traced by galaxy clusters to redshifts 
z = 1 and beyond. The XMM-LSS sensitivity limit is ~ 1000 
times deeper than REFLEX i.e. ~ 4 X 10~^^ ergs~^cm~^ in the 
[0.5 - 2] keV band for point sources. Moreover, the XMM-LSS is 
able to make a systematic exploration for massive clusters out to at 
least z ~ 1.5 |Refregier et1l]p002| ), [Bremer et al.| |2006 ), Pa caud] 
|et al.| ( [2006] )]. Another key improvement is the spatial resolution of 
XMM: ~ 6" FWHM on axis compared to ~ 20'' for ROSAT. 

The two major requirements of the X-ray processing were to 
reach the sensitivity limit of the data in a statistically tractable man- 
ner in terms of cluster detection efficiency, and so to provide the 
selection function of the detected objects. To achieve these goals, it 
was necessary to design a new two-step X-ray pipeline, combining 
wavelet multi-resolution analysis with maximum likelihood fitting; 
Poisson statistics being used in both steps ( jPacaud et al.|2006 ). We 
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stress here that at our sensitivity and spatial resolution, the survey 
is primarily limited by surface brightness, rather than by flux as 
assumed by past generations of surveys. This led us to define clus- 
ter selection criteria in a two-dimensional parameter space, corre- 
sponding to specific levels of contamination and completeness, as 
discussed below. 



Building on the preliminary results from the first square de- 
gree of the survey ( [Pierre et al. !2Q06 ), we have defined a well- 
controlled cluster sample covering all currently available XMM- 
LSS observations, i.e. about 5 deg^. The sample comprises some 
30 clusters with fluxes in the [0.5-2] keV band ranging from 1 to 
50 xlO"^"^ erg s~^cm~^. They are referred to as 'First Class', or CI 
clusters, because the criteria used to construct the sample from a 
two-dimensional X-ray parameter space guarantee no contamina- 
tion by point-like sources. The observations were performed in a 
rather homogeneous way (10-20 ks exposures) and enable us, in 
addition, to estimate temperatures for all the sample clusters. This 
is a key datum which, along with careful modelling of the survey 
selection efl'ects, allow us to tackle for the first time on real data 
the calibration of the mass-observable relations in a cluster survey 
- something which has been mostly addressed in a formal way thus 
far [e.g. [Majumdar & M ohr' ('2003') and references therein]. This 
leads us to quantitatively investigate some of the degeneracies af- 
fecting the cosmological interpretation of the redshift distribution, 
dn/dz, of cluster number counts. 




37.0 36.5 36.0 35.5 35.0 

Figure 1. Sky distribution of the CI clusters. The colour scale indicates 
the on-axis exposure time of each pointing in seconds (after particle flare 
filtering). Pointings marked by a cross are strongly affected by flares and 
will be re-observed during AOS. 



2 CONSTRUCTION OF THE SAMPLE 



The paper is organized as follows. In section 2, we summarise 
the principles of our X-ray analysis, provide a detailed calculation 
of the selection function, and subsequently present the cluster sam- 
ple. The next section describes the spectroscopic conflrmation of 
the clusters, and the determination of their X-ray temperature and 
luminosity. Section 4 provides a thorough discussion of the Lx -T 
relation: special care is given to the modelling of the selection ef- 
fects in the analysis of possible evolutionary trends. This is an espe- 
cially informative exercise as it is the flrst time that the behaviour of 
the relation has been explored for a well controlled sample. Section 
5 is devoted to the modelling of the sample; starting from the linear 
spectrum of the initial density fluctuations, we compute the dark 
matter halo mass spectrum and infer the halo X-ray luminosity and 
temperature using empirical scaling laws compatible with our data. 
The cluster population thus obtained is then folded through the 
XMM-LSS selection function. The derived dn/dz is compared with 
that observed, and the model provides us with mass-observable re- 
lations. The model also allows us to further explore (Section 6) 
the impact of several cluster parameters on classical cosmological 
tests; for instance, we examine the sensitivity of our observed dn/dz 
to the shape of the M-T relation in the group regime, to the amount 
of scatter in the scaling laws, and to changes in cosmological evo- 
lution resulting from diff'erent equations of state for dark energy. 



2.1 The dataset 

To date, the X-ray data received consist of a single region of 
roughly 6 deg^ covered by 51 XMM-Newton pointing^ Most of 
the observations (32 pointings, hereafter B pointings) were ob- 
tained through XMM A01/A02, and have a nominal exposure time 
of 10^5, whilst 19 (hereafter G pointings) are deeper 2 x lO'^^ guar- 
anteed time observations. The latter constitute the XMM Medium 
Deep Survey (XMDS, piiappetti et al.|2005l l. 

Calibrated event lists were created using the XMM- SA^ tasks 
emchain and epchain, and periods of intense proton flares were 
filtered out using the method suggested by [Pratt & Arnaud|p002| l. 
Out of the 51 observations, 6 (including one G pointing) were too 
strongly contaminated to properly monitor and remove the high 
flare periods. These pointings are scheduled for re-observation dur- 
ing XMM A05 and, consequently, are not included in the analysis 
of the present paper, which covers only 5 deg^. 

Details of the observations can be found in the p ublication of 
the full X-ray catalogue pertaining to this first 5 deg^ ( Pierre et al. 
|submitte"d| ). 



2.2 Data pre-analysis 

For each pointing, images were generated in the [0.5-2] keV band 
from the filtered event lists using the XMM-SAS task evselect. 
These were subsequently analyzed using the XMM-LSS 



Throughout this article, data are analyzed using the cosmo- 
logical parameters estimated by [Spergel et al.| ( |2006| ), namely: 
Ho = 73 km s-^ Mpc'K = 0.24, Qa = 0.76, = 0.041, 
n, = 0.95 and os = 0.74. 



^ These are currently being complemented during XMM-Newton AOS so 
as to cover the full 10 deg^ of the Spitzer/SWIRE area surveyed in the 
XMM-LSS field. 



XMM Science Analysis System, http://xmm.vilspa.esa.es/sas/ 
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source detection pipeline (described in detail by |Pacaud et al.| 
[2006 ): First, the images are filtered using a wavelet multi-resolution 
algorithm ( Starck & Pierre, 1998 ) that was specifically designed to 
properly account for the Poisson noise in order to smooth the back- 
ground and lower the noise level, while keeping unchanged the rel- 
evant information. Then SExtractor is used on the filtered images 
to extract a very deep primary source catalogue. Finally, detailed 
properties of each detected source are assessed by Xamin, a maxi- 
mum likelihood profile fitting algorithm that we have developed to 
characterize extended sources in XMM-Newton images. 

2.3 Source selection process 

2.3.1 Definition of the sample 

As shown by |Pacaud et al.| { [2006 1, the probability of detecting an 
extended source at the XMM-LSS sensitivity and resolution de- 
pends on both its flux and angular extent. This led us to abandon 
the simple concept of a flux limited cluster survey and to define 
a system of extended source classes: the final source selection is 
performed in the Xamin output parameter space, and several sam- 
ples are defined, allowing for various amounts of contamination 
from point sources and spurious detections. The sub-samples, or 
classes, have been defined by means of extensive 'in situ' simula- 
tions, involving a large range of cluster fluxes and apparent sizes, 
for the nominal exposure of the XMM-LSS, i.e. 10^5. In this paper, 
we focus on the Class 1 sample (hereafter CI) which was defined 
from our simulations to be the largest uncontaminated extended 
source sample available. It is obtained by selecting candidates with 
extension > 5'', extension likelihood > 33, and detection 
likelihood > 32. Given that no spurious extended detection was 
found in our whole simulated dataset, we are confident that the 
false detection rate is close to zero, and visual inspection of the 
X-ray/optical overlays for the CI sources showed the existence of 
a cluster of galaxies or of a nearby galaxy as an optical counterpart 
in every case. Furthermore, we find that the CI sample coincides 
with clusters for which a reliable temperature can be obtained us- 
ing X-ray spectroscopy from the survey data (section [T2| . 

Detailed information on the CI selection process can be found 
in [Pacaud et al.| { [2006| ). 

2.3.2 Selection function 

Given the varying eff'ective exposure time per pointing, the survey 
coverage is not uniform (7 < texpl^^^s < 20, Fig.[lJ. Consequently, 
similar likelihood values do not exactly pertain to the same objects 
across the survey. Rather than lowering all exposures to 7 x lO^s, 
we apply the CI criteria to every pointing in order to maximize the 
size of the cluster sample. The probability of detecting a cluster of 
a given flux and extent for any pointing is subsequently derived by 
applying an analytic correction to the \(fs simulations, scaling as 
a function of exposure time the S/N produced by such an objecj^ 
This is justified by the fact that the range spanned by the exposure 
times is modest (1/2 < texpl^^'^s < 2), implying also that source 
confusion does not change significantly. The survey selection func- 
tion, i.e. the probability to detect a source of a given flux and appar- 
ent size as a CI cluster, is obtained by integrating the contributions 
from all pointings. It is this function that is used in the following 



^ The mean background values measured by [Read & Ponman] j2003) are 
assumed 




0.02 0.04 0.06 0.08 0,10 

M1+M2 + PN count-rate (cts/s) 

Figure 2. The survey coverage is displayed in a 2D parameter space: the sky 
coverage is a function of both cluster flux and extent. The dashed lines are 
the result of extensive lOks simulations. The slightly shifted white lines are 
the analytical corrections accounting for exposure variations across the sur- 
veyed area, hence indicating the effective selection function of the current 
data set. The sample completeness corresponding to each line is indicated 
in terms of both the sky coverage and the percentage of detected sources. 
Extent values correspond to the core radius of a yS-model with j3 = 2/3. The 
count-rate to flux conversion assumes a 2 keV spectrum at z = 0. 

for cosmological applications. It is displayed in the form of a sky 
coverage in Fig. [5] The maximum area covered reaches 5.2 deg^. 
This happens e.g. for sources with 20''core radius, and fluxes above 
5 X 10"^"^ erg s~^cm~^. Our sensitivity drops by roughly a factor of 
two for sources with flux around 1.5 x 10"^"^ erg s~^cm~^ and 20'' 
core radius, or 5 x 10"^"^ erg s~^cm~^ flux and 75'' core radius. 

2.4 The sample 

Applying the CI selection criteria to the inner 11. 5 'of the 45 valid 
observations yields 37 detections, among which 4 are duplicated 
sources detected on adjacent pointings. Optical/X-ray overlays re- 
veal that 5 of the remaining 33 sources are actually nearby X-ray 
emitting galaxies whose properties are presented in Appendix |C] 
Hence we have a total of 28 galaxy clusters. An additional unam- 
biguous CI source is detected in pointing B13 (highly contami- 
nated by flares); we provide the basic properties of this source in 
the present article, but it is not used in the scientific analysis as we 
cannot accurately compute the selection function for this pointing. 
The cluster sample is listed on Table [T] We note that a complemen- 
tary cluster sample of about the same size has also been identified 
at lower significance in the same area; these clusters belong to our 
C2 and C3 classes, for which the selection criteria are much less 
well defined. They will be published in a separate paper (Adami et 
al., in prep.). 



3 ASSESSING INDIVIDUAL SOURCE PROPERTIES 
3.1 Spectroscopic validation 

Each CI X-ray source was confirmed as a bona-fide galaxy cluster 
at specific redshift via optical photometry and spectroscopy. 
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The overall procedure used to confirm individual clusters is 
very similar to the approach taken in previously published XMM- 
LSS p apers (Valtchanov et al. 2004l|Wims et al.|2005, Pierre et al.| 
|2006| ) and is summarised below. 

A c ombination of either CTIO/MOSAIC II Rz' ( |Andreon et| 
al.|2004 | or CFHT/MEGACA]V0wgr/z imaging was used to asso- 



ciate the location of each X-ray source with the spatial barycentre 
of a significant overdensity of galaxies displaying characteristically 
red colours. Based on the CTIO data, galaxies lying within a given 
colour tolerance of this "red sequence" were then flagged as candi- 
date cluster members and given a high priority in subsequent multi- 
object spectroscopic observations. A small number of low X-ray 
temperature (< 1 keV) groups at moderate redshift (z > 0.2) and 
local {z < 0.2) compact groups did not display a statistically signifi- 
cant red galaxy overdensity. These systems were inspected visually 
and spectroscopic targets were assigned manually. 

Due to the moderately large size and the extended redshift 
range covered by the sample, cluster targets were observed using 
a number of facilities over several observing semesters. Details of 
the observing configurations can be found in Table |2] 

Spectroscopic data were reduced using standard procedures 
described in detail in previously published XMM-LSS publica- 
tions ( Valtchanov et al.|2004||Winis et al.|2005l [Pierre et al.|2006| ). 
Where possible, redshift values and associated uncertainties of in- 
dividual galaxy targets were computed via cross-correlation with 
galaxy reference templates. In the remaining cases redshifts were 
assigned manually, from emission features. 

The nominal condition adopted by XMM-LSS to confirm a 
cluster redshift is to observe three concordant redshifts (typically 
Az < 0.01) within a projected scale of about 500 kpc of the cluster 
X-ray centre. Only cluster XLSSC-052 does not fulfill this condi- 
tion: the X-ray emission is associated with a pair of z = 0.06 galax- 
ies very close to a bright star. Except for this specific system, the 
number of galaxy members confirmed per cluster ranges between 
3 and 15, with a typical cross-correlation velocity error of order 
50-150 km s~^ per galaxy. Despite this velocity accuracy, we are 
limited for a large fraction of our sources by small-number statis- 
tics and consequently quote redshifts only to +0.01. In each case, 
the final cluster redshift is computed from the unweighted mean of 
all galaxies lying within +3,000 km.s~^ of the visually assigned 
redshift peak. 



3.2 Spectral analysis 

In order to measure the temperature of the intra-cluster gas. X-ray 
spectra were extracted in a circular aperture around each source. 
The corresponding background emission was estimated within a 
surrounding annulus having inner radius large enough for the clus- 
ter contribution to be considered negligible. Preliminary modelling 
of the cluster surface brightness profile allowed the determination 
of the optimal extraction radii in terms of S/N. 

The resulting spectra were fitted using to a thermal 

plasma model (APEC) assuming a fixed hydrogen column density 
set to the Galactic value as derived from H I observations by Dickey" 
|& Loc km"anl (l990 ). The metal abundance of the gas was held fixed 
during the fitting process at 0.3 times the solar abundance, as esti- 
mated by |Grevesse & Sauval| ( |1998) . As explained in [Willis et al.| 



( |2005| ), the cluster spectra were constructed imposing a minimum 
requirement of 5 background photons per bin in order to avoid the 
apparent bias we identified in Xspec temperature estimates when 
using the Cash statistic on very sparse spectra. Our simulations 
( [Willis et al.|2005 ^ showed that this procedure provides quite reli- 
able temperature measurements (+10 - 20%) for ~l-3 keV clusters 
having only a few hundred counts. We further investigated the im- 
pact of fixing the metal abundance at 0.3 Z©, by computing best fit 
temperatures obtained using extreme mean abundances of 0. 1 and 
0.6 Zq. In most cases, the temperatures fell within the Icr error bars 
from our initial fit. For five systems (namely XLSSC-008, XLSSC- 
028, XLSSC-041, XLSSC-044 and XLSSC-051) one of these two 
extreme measurements just fell a few percent outside of our error 
range. The measured temperatures are presented in Table [T] 



3.3 Spatial analysis 

3.3.1 Surface brightness modelling 

To accurately determine the cluster fluxes and luminosities we 
modelled the observed photon spatial distribution in the [0.5-2] keV 
band with a radial /^-profile: 



S{r) 



(1) 



[l+(r//?,)2]3^-i/2' 

using a refined version of the method described in [Pierre et"aL] 
([2006}. 

For each source, we started by fitting the mean background 
levels (vignetted and particle components) over the inner 13' of the 
pointing, excluding all sources detected by the pipeline. Then, for 
each EPIC camera, the photons were binned within y annuli cen- 
tred on the cluster peak. The resulting profiles were subsequently 
rebinned imposing a minimum S /N of 3 with respect to the esti- 
mated background level, weighted by the annular exposure times, 
and finally co-added to build a M0S1+M0S2+PN count-rate pro- 
file. The fitted model is constructed by convolving the circular /3- 
profile w ith an analytical par ametrization of the P SF ( [Ghizzardi et[ 
al.|200l t, as implemented by [Arnaud et al.[ ( [2002| l. The;^^^ statistic 
is computed over a discrete grid of j3 and Rc values, with the value 
of the normalization coefficient, Sq, optimised analytically. 

As already discussed in [Pierre et al.[ ( [200 6 ), the majority of 
the clusters are faint, and the detected photon distribution in many 
cases represents only a fraction of the extended X-ray surface 
brightness distribution. Under such conditions the parameters j3 and 
Rc are degenerate when fitted simultaneously, limiting the extent to 
which 'best-fitting' parameters can be viewed as a physically re- 
alistic measure of the cluster properties, although they provide a 
useful ad hoc parametrisation (Appendix [A|. For this reason, we 
do not quote here best fitting values of /3 and Rc derived for each 
confirmed cluster. A dedicated analysis of the mean cluster profiles, 
obtained by stacking the data as a function of redshift and temper- 
ature, is underway (Alshino et al in prep.). 

The photon count-rate within a specified radius is obtained by 
integrating the best-fitting spatial profile. The conversion into [0.5- 
2] keV flux and unabsorbed rest-frame bolometric ([0.001-50] keV) 
luminosity is performed via Xspec using the cluster temperatures 
previously derived. 



^ Data are taken from the CFHT Wide Synoptic Legacy Survey. See the 
URL www. cfht .hawaii . edu/Science/CFHTLS/ for further details. 



http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 



3.3.2 Luminosity and flux determination 

Luminosities are integrated within rsoo i.e. the radius at which the 
cluster mass density reaches 500 times the critical density of the 
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Table 2. Instrumental characteristics for each spectrograph configuration employed during the observations. All spectral 
observations were performed with a slit width of between 170 and 174. 



Telescope Instrument Grism + Filter Wavelength Pixel sampling Spectral Identifier 

interval (A) (A pix~ ^ ) resolutioip] (A) 



VLT 


FORS2 


300V + GG435 


4000-9000 


3.2 


14 


1 


VLT 


FORS2 


600RI + GG435 


5000-8500 


1.6 


7 


2 


VLT 


FORS2 


600z + OG590 


7500-10000 


1.6 


7 


3 


Magellan 


LDSS2 


medium-red 


4000-9000 


5.1 


14 


4 


NTT 


EMMI 


Grism #3 


4000-9000 


3.0 


8 


5 



Estimated via the mean full- width at half-maximum of the bright, isolated arc emission lines. 



Universe at the cluster redshift. As in ' Willis et al.| ('2005^ and'Pierre] 
|et al.| ( [2006| , this radius is estimated from the cluster mean temper- 
ature using the M-T relation of |Finoguenov, Reiprich & Bohringer| 
( |2001| ), converted to AC DM cosmology, which gives: 



rsoo = 0.375 T^-^^ h.^izT' Mpc, 



(2) 



where T is expressed in keV and hj2 is the bubble constant in units 
73 km s~^Mpc~^ Although "total" fluxes are often quoted in the 
cluster literature, our simulations show that the present data do not 
allow us to reliably perform such measurements (see Appendix[A|. 
To limit extrapolation uncertainties, aperture flux values are com- 
puted by integrating within a fixed radius of 0.5 Mpc. As conspic- 
uous in Table [T] 0.5 Mpc is generally smaller than 2 x Rfu except 
for the nearby low temperature groups XLSSC-039, XLSSC-021 
and XLSSC-052. We chose to compute the flux within 0.5 Mpc 
(rather than Rsoo, which is generally smaller), since this is a simi- 
lar approach to aperture photometry. For all clusters it corresponds 
to about 2/3 of the total flux, i.e. integrated to infinity, assuming 
a profile defined by /5 = 2/3, Rc = 180 kpc. We emphasise that 
these fluxes are only used in the comparative analysis of the log(N)- 
log(S) relation, while L500 and M500 (along with T) are the actual 
physical quantities used in our cosmological modelling and in the 
subsequent discussion. 

Finally, the Icr errors on the extrapolated fluxes and luminosi- 
ties are computed by identifying the region of the (Sq, Rcfi) param- 
eter space where ^ 1^ computing the extreme values of the 
extrapolated count-rates allowed by these models (see Appendix [A] 
for further discussion on the flux measurement accuracy). 



3.33 Mass determination 

Based on our best spatial fit profile, we estimate cluster masses, as- 
suming that the gas is isothermal (the limited number of photons 
does not allow us to derive temperature profiles). Under these as- 
sumptions, for M500 in Mq, the hydrostatic equilibrium assumption 
yields (see e.g. |Ettori|2000l ): 



M. 



1.11 XO^^xpRJ 



^ ^ -^500 



(3) 



where X500 = r^oo/Rc (with r^oo given by Eq.|2]), R^ is expressed in 
Mpc, and T in keV. 

Values of rsoo, flux and luminosity for the CI clusters are listed 
in Table [1] together with M500. 




1.2 



Figure 3. Redshift distribution of the CI sample. The colour scale indicates 
the cluster mean temperature for each bin (unweighed mean of the individ- 
ual cluster temperatures in keV) . 



4 RESULTS 

4.1 Global properties of the sample 

The detected clusters span the [0.05-1 .05] redshift range with a pro- 
nounced peak around z ~ 0.3 (Fig. [3]). Roughly half of the objects 
have a temperature lower than 2 keV, pointing to a large fraction of 
groups in our sample. As a natural consequence of its sensitivity, 
good PSF and dedicated source detection procedure, the XMM- 
LSS systematically unveils for the first time the z ~ 0.3, T < 2 keV 
cluster population on large scales. We note also that none of the ob- 
jects detected so far exhibits strong lensing features. In the follow- 
ing, we will generally use the word "cluster" regardless of the tem- 
perature, while "group" specifically refers to objects below 2 keV. 
A noticeable feature of the sample, illustrated by Fig. [3] is the in- 
creasing mean cluster temperature with redshift. Given the strong 
correlation between X-ray luminosity and temperature (Sec. |4.3| ) 
this can readily be understood as the unavoidable Malmquist bias. 
Moreover, the current sky coverage of the survey makes the sample 
likely to be aflfected by cosmic variance. Both issues will be given 
special attention in the following cosmological analysis. 

Additional information on individual sources, including possi- 
ble contamination by AGN emission, are discussed in Appendix [B] 
We also display the optical images of each cluster with confirmed 
member galaxies and X-ray contours overlayed. 
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Table 3. Trends in the literature for the modeUing of F(z) as a power 
law of (1 + For comparison, constraints from the CI sample were 
added, ignoring the selection effects. 



Ref. 



power of (1+z) 



1 2002 



Vikhlinin et al. ( 2002| 

Nm ^icki, Sornig, & Henry 
Ettb ri et al. ( 20^| 
Lumb et al. ( 2QCT 
Kotov & Vikhlinin ( 2005^ 
Maughan et al. ( 2006 ) 
CTclusters with 0.1 < z < 0.4 (14 sources) 
CI clusters with 0.4<z<l.l(10 sources) 
All CI clusters above z = 0.1 (24 sources) 



1.5 ±0.3 
2.1 ± 1.1 
0.6 ± 0.3 
1.5 ±0.3 
1.8 ±0.3 
0.8 ± 0.4 
2.3 ±0.8 
1.3 ±0.5 
1.5 ±0.4 



4.2 The Msoo-Lx relation 

The correlation between cluster luminosities and masses derived 
from the spatial analysis is shown in Fig.|4] Once self-similar evo- 
lution is assumed, our data points appear to be continuous with the 
massive cluster samples of [Zhang et al.| ( [2006t and [Zhang et al.| 
( [200 7 ) [respectively at z ~ 0.3 and z ~ 0.2]. Moreover, the scatter 
in the relation is relatively low (i.e. almost comparable to that of 
the high mass sample). This overall consistency is remarkable in 
the sense that, for half of the clusters, less than 500/300 photons 
were available for the spatial/spectral fitting. This provides an im- 
portant indicator of the reliability of our measurements, and could 
be further interpreted as adding weight to the self- similar evolution 
hypothesis. 

The comparison between our data points and the indicative 
M-Lx relation derived from [Arnaud, Pointecouteau & Pratt[ ( |2005| ) 
and [Amaud & Evr^d| ( |1999| ) also suggests some flattening of the 
correlation at the low mass end. Our luminosity sampling, how- 
ever, suffers significant Malmquist bias, and masses derived from a 
single isothermal /?-model fitting are known to be underestimated 
( [Vikhlinin et al|2006| . Both effects tend to contribute to the appar- 
ent flattening in the M-Lx relation. 



4.3 The Lx-T relation 

The correlation between the observed luminosity and the tempera- 
ture of our clusters is show n in Fig. [5] The co mparison with the 
local relation measured by [Arnaud & Evrard[ fl999 ) suggests a 
positive redshift evolution of the luminosity at a given tempera- 
ture. To test this, we performed an analysis of the enhancement 
factor, F(z, T) - Lx/Lx(T,z = 0), with respect to this local refer- 
ence. In most of the intra-cluster gas models, F is a simple func- 
tion of redshift: the self-similar assumption yields F(z) = E(z) 
where E(z) = H{z)IHq, is the evolution of the Hubble constant. 
More elaborate models that include non-gravitational physics (i.e. 
cooling and heating) generally propose an evolution of the form 
F{z) = (1 + zT (e.g. |Voit|200 5). All these models assume that the 
entire cluster population follows a unique evolutionary track. In the 
present study, we adopted this hypothesis, since the Malmquist bias 
prevents us from investigating the extent to which evolution could 
be a function of cluster temperature. 

To date, various attempts to constrain F{z) as a power law in 
(1 + z) have yielded discordant results (see Table [3]). Several ex- 
planations for this discrepancy have been invoked: poor statistics, 
deviations between the several available local reference relations, 
disparity in mass and redshift of the current high-z cluster samples 
(if F is not a simple function of z), as well as biased samples. Al- 
though limited by photon statistics, the XMM-LSS cluster sample, 
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This work, locol compoct groups 
This work, contominoted sources 
This work, other sources 



• •• 



* 

* * 




xE(z) [h-^ M, 



Figure 4. The Mass-Luminosity relation for the sample clusters. Because 
of the large redshift range spanned by the data, the mass and luminosity 
parameters are scaled assuming self-similar evolution (factor Eiz), see Sec. 
[4.3) . For comparison, the massive cluster samples from |Zhang et al.|{T006) 
and 'Zh ang ern?| ( |2007) are also shown. The dashed and solid lines also 
show the mean Msqq-Lx relation infered from the M^qq-T and Lx-T mea- 
sured in [ Arnaud, Pointecouteau & Pratt|p005i and | Arnaud & Evrard|(1999) 
respectively for T > 3.5 keV and T > 2 keV 



with its well controlled selection function, provides an important 
opportunity to test the impact of the selection process in such stud- 
ies. 



In a preliminary analysis,'Willis et ar] ( [2005 l measured a mean 
value F - 1.46 from an initial XMM-LSS sample covering the 
[0.30 - 0.43] redshift range. Another estimate using the Dl sub- 
sample (.Pierre et al. 2006 section 3.4) also found a luminosity en- 
hancement for the z ~ 0.3 domain, and showed that it could be 
statistically significant only above T ~ 1.5 keV. However, none of 
these studies fully addressed the impact of selection biases on the 
derived results. 

In the present article, involving a larger sample, we perform a 
self-consistent likelihood analysis of our data parameter space, for 
a given F{z). We excluded from the analysis XLSSC-050 (contami- 
nated by flares), XLSSC-018 (contaminated by an AGN), and the 3 
local groups (z < 0.1 and T < \ keV) for which our measurements 
could be affected by a large radial extrapolation and by X-ray emis- 
sion from member galaxies. The available sample for this analysis 
is thus restricted to 24 objects. 



4. 3. 1 The likelihood model 

Folding our selection function into the analysis of the Lx-T re- 
lation cannot be achieved by the usual fitting. Instead, we have 
to estimate the likelihood of our data, given a F{z) model. We thus 
start by assuming that cluster luminosities are distributed by a log- 
normal probability density function (PDF) around the mean value 
Lmean(T,z) = L(T,z = 0) X F(z), which we will call the true distri- 
bution: 



dPtrue(Lx\T,z) 

dlnLv 



1 



^/27TC 



■ exp 



\nLx\T 



In (Lx/LfneanY 
^^\nLx\T 



(4) 
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T (keV) 



Figure 5. The L500-T relation. Same symbols as in Fig.|4] The sources that 
are used for the analysis of the Lx - T evolution (circles) have redshifts in 
the range [0.14-1.05] and fall on average above the local Lx-T relation of 
[Arnaud & Evrard|jl999t (shown by the solid line). 

An estimate of the scatter in the local Lx-T relation is provided by 
[Stanek et al.|(|2006[). Th e authors fitted a single dispersion value 
to the Reiprich & Bohringer^ (2002 ) sample (thus assuming that it 
is the same from the group regime to large clusters) and obtained 
cTinTiLx - 0-25. We further assume that this scatter does not evolve 
with redshift. Combining the lognormal distribution with the power 
law shape of the mean Lx-T relation, we get crinLj^ir ~ 2.7x0.25 ~ 
0.7 where 2.7 is the approximate Lx-T slope in the [Reiprich &| 
|Bohringer| ( [2002l ) sample. 

The normalised PDF for detection of a cluster with such Lx 
and T is then the product of the 'true' PDF with our selection func- 
tion: 

dPdet{Lx\T,z) 

dlnLx 
where: 



PoiT,z)x 



dP,„e(Lx\T,Z) 

dlnLv 



xf(Lx,T,z), 



(5) 



(I 



dP,rue(Lx\T,z) 



xf(Lx,T,z) 



dilnLx) 



(6) 



dlnLx 

We compute f(Lx, T,z) (from Fig [2]), assuming a canonical value 
of 180 kpc for the cluster core radii. The resulting modification 
of the PDF is illustrated in Fig. [6] at diff'erent redshifts assuming 
self- similar evolution: including selection eff'ects in the likelihood 
model renders underluminous cool clusters undetectable. Note that 
for clusters with temperatures close to the detection threshold at 
any given redshift, the bias in mean Lx at a given T due to selection 
eff'ects can be a factor 2 or more. 

Practically, errors on the luminosity can be considered as neg- 
ligible, since the fractional errors on temperature are usually much 
larger and the Lx-T relation goes like T^ . The likelihood of an 
observed cluster C being drawn from a given F{z) is thus obtained 
by marginalizing over the possible cluster temperatures, i.e.: 



L{F\C) : 



dPi„(Lx\T,z) 



Jin Ly 



X Pc{T)dT, 



(7) 



where Pc{T), the temperature PDF of the cluster, is derived from 
the C- statistic distribution provided by the Xspec's steppar com- 
mand. Finally, the likelihood of the enhancement factor is com- 




0.5 as D.7 o.g o.g i,d 



4.0 5.0 &.0 7.0 



T (keV) 

Figure 6. Comparison of the modelled Lx-T distribution with and without 
selection effects for several redshifts, assuming self-similar evolution. The 
light and dark shaded regions show, for each temperature, the luminosity 
interval than contains respectively 95.5 and 68.3% of the expected detec- 
tions. The solid line gives the maximum of the luminosity PDF for each 
temperature; it is cut at the temperature for which the detection probability 
(integrated over Lx) falls to 5%. For comparison, the local Arnaud & Evrardj 



jl999) relation, evolved according to the self-similar model, is indicated by 
the short dash line, and the 1 and 2cr bounds arising from our scatter model 
are shown by dotted lines. The long dash line is the local relation. 



puted as the product of the detection probabilities, for the given 
F(z), of all clusters pertaining to our sample: 



L(F) = Y\l(F\Q) 



(8) 



To estimate errors on fitted parameters, we identify -2 InL(F) 
with SLX^ distribution, as is asymptotically valid in the limit of large 
samples according to Wilks' theorem, and quote 68% confidence 
intervals. The same identification is used to assess how much a 
given model deviates from the best fit F(z) by correcting the mea- 
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Figure 7. Evolution of the Lx-T relation, (a): Raw analysis. The data points from several studies are plotted (stacked into redshift bins for clarity). Whenever 
possible, we converted luminosity estimates from the other samples to L500 using the information provided by the authors. The following differences remain: 
the luminosities of th e|Lumb et al.|{2004) clusters are estimated within the virial radius; the data points of |Maughan et al.|{2006) also include the clusters 
from |Vikhlinin et aL]j2002} and have luminosities estimated within R200. Overlaid are several enhancement factor fits from our baseline analysis: the (1+z) 
power law and ad-hoc 2 parameter fits to our data alone (resp. short and long dashed lines), and the (l+z) power law and ad-hoc 2 parameter fits (resp. 
dot-dashed and 3dot-dashed lines) fitted to the CI clusters combined with the samples of Kotov & V ikhlin'iii|j20 05'> and Ettori et al. (2004 1. For comparison, 
the self-similar prediction is indicated by the solid line, (b): Taking into account selection effects. The filled circles and short dashed line recalls the location 
of our raw data points, and best fit model from the preliminary analysis. The corrected enhancement factor fit for the 2 parameter model is shown as the long 
dashed line. The final 1 parameter fit and its Icr confidence interval is displayed by the solid line and the shaded region. Expectations from several intra-cluster 
gas models are plotted for comparison: the self-similar predictions (dotted line) and two models by |Voit|{2005) including non gravitational physics (dot-dashed 
and 3dot-dashed lines) 



sured I^x^ to the equivalent value that yields the same probability 
for 1 free parameter. 



4.3.2 Rawfit 

In a first step, we neglect the selection function in our formalism 
(i.e. effectively imposing f{Lx, T,z) - I in the likelihood compu- 
tation). This will both provide us with a reference point to assess 
the impact of the selection process, and allow for comparison with 
previous work, where the issue of sample selection has been ig- 
nored. 

As shown in Table |3] fitting the usual F{z) - (1 + z)'^ yields 
a = 1.5 + 0.4, which is consistent with [Vikhlinin et al.| ( |2002| ), 
ILumb et al.| ( [20041 ) and |Kotov & Vikhlinin| (^200?). Nevertheless, 
this simple model does not seem to correctly reproduce the obser- 
vations over the full redshift range: the values of a obtained over 
[z = 0.1 - 0.4] and [z = 0.4 - 1.1] are incompatible at the one 
sigma level (see Table |3]). Interestingly, [Ettori et al.| ( [2004| noticed 
a similar discrepancy within their own sample, but did not consider 
its implications for their subsequent analysis. 

In Fig. [7^, we display simultaneously our data points along 
with those from previous studies. Error bars are large, but the 
points suggest a non-monotonic evolution, with a maximum around 
z ~ 0.5, thus excluding an enhancement factor of the form F(z) = 
(1 + zT- This may explain the diverging results of Table |3] Con- 
sequently, we fitted an ad hoc two-parameter model of the form 
F(z) = (1 + zT X E(z)^. This model is intrinsically degenerate and 
would need very accurate data to individually constrain a and but 
has the ability to reproduce the apparent non-monotonic evolution. 
Since cases a - - 1, and p - correspond to the self-similar 



and power law models respectively, we may evaluate through this 
parametrization how much the observations deviate from any of 
them. 

A fit over the CI clusters yields a - 4.7 and p - -5.4 and 
excludes the self-similar evolution at the 3.5cr level. Our best (1 + 
zT model is however less than 0.4cr away from the best fit and the 
evidence of a non-monotonic evolution is thus rather weak from our 
data alone (the probability of such a deviation being real is roughly 
45%). This is not the case any longer wh en in cluding also in the fit 
the data from |Kotov & Vikhlinin| ( |2005| ) and [Ettori etah] ( |2004|' ^ " 



for which we could convert the luminosity to L500. A simple power 
law model is then ruled out at the 8.4cr level based on 59 sources. 



4.3.3 Taking into account selection effects 

With this comparison baseline, we can now turn to investigate the 
impact of our selection process on the observed Lx enhancement 
(restoring the correct value of f(Lx,T,z) in the likelihood). This 
task is tractable here because the XMM-LSS CI selection function 
has been thoroughly assessed. 

Maximizing the likelihood over our ad hoc two-parameter 
model now yields or = 1.3 and p = -1.3. As can be seen in 
Fig. [TJd, the inferred evolution rate is considerably lower than with 
the raw fit. As a result, the self-similar expectation is now less than 
7 X 10~^cr away from the best fit model, and the statistical evi- 
dence for a non-monotonic evolution is virtually null. To evaluate 
the remaining deviation from self- similarity, we switch back to a 

^ We assumed a gaussian PDF of the cluster temperature for these sources 
in our likelihood model 
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function from |Bardeen et al.] ( |1986| ), taking into account the shape 
parameter: 



Figure 8. Similar to figure|6]for z = 0.3, with the 0.2 < z < 0.4, with CI 
clusters overlaid. Again, the dashed line is the evolved model oflArnaudX 



Evra rd (1999 1. The triangles indicate the Dl clusters presented by Pierre et 
arn2006 >. The only cluster (XLSSC-018, square) that falls well below the 
Arnaud & Evrard (1999 1 Lx-T relation, and at the limit of our 2cr contour. 



is likely to be contaminated by an AGN (see Appendix [b). 



one parameter model by fixing p - \. The best fit value for a is 
then -0.07^Q 53, and the corresponding F{z) is shown in Fig.jrjD. 

This latter analysis points to a negative evolution of the L^-T 
relation (i.e. L{z\T) increases with z) over the < z < 1 range. 
These results quantitatively favour the self-similar hypothesis, al- 
though our current data set is still marginally consistent with no 
evolution at all as is clear from Fig.|7]3. 

This is illustrated in more detail in Fig. |8] for the z ~ 0.3 in- 
termediate redshift, where the 1-2 keV groups make up the peak of 
the XMM-LSS sensitivity (cf. sectionjS]). 

These results are still preliminary because of the poor statis- 
tics (in both the number of sources and the temperature estimates). 
Moreover, they are likely to be very sensitive to the choice of 
the local Lx-T reference (including scatter) due to the redshift- 
temperature correlation within the sample. Nevertheless, our find- 
ings seriously question any attempt to assess the evolution of clus- 
ter scaling laws without modelling the full source selection process. 

In the rest of the article, we assume that the cluster scaling 
laws evolve following the self- similar model, as suggested by the 
above analysis. 



5 COSMOLOGICAL MODELLING 

In this section, we compare the observed properties of the CI sam- 
ple with cosmological expectations based on the latest WMAP con- 
straints ( |Spergel et aL|2006| ). 



5.1 The model 



We model the cluster population and its evolution following' Pacaud| 
|et al.| ( [2006] ), but with slightly modified assumptions as to the clus- 
ter scaling laws. 

The linear power spectrum P{k) of dark matter fluctuations is 
computed using an initial power law of index k^' and the transfer 



r = ^l^h X exp [-at (l + V2/z/n^)] 



(9) 



introduced by |Sugiyama| ( |1995t . The overall power spectrum is nor- 
malised at z=0 to cTg, and the redshift evolution is obtained from the 
linear growth factor ( [Carroll, Press, & Turner 1992 ). Based on this 
power spectrum, we use the |jenkins et al., ( ,200 1 ) formalism to de- 
rive the dark matter halo mass function dn/dm(z). This procedure 
is an accurate fitting of the mass function obtained from numerical 
simulations, provided one defines the mass of the halo to be that in- 
cluded inside r2oob^ the radius enclosing an overdensity of 200 with 
respect to the mean background density. 

In order to reproduce the CI selection function, we need to 
translate the mass parameter into observable parameters, which is 
practically achieved by means of scaling laws. Usually, halo masses 
in such relations are defined with respect to the critical density of 
the Universe (and not to the background density). Hence, the need 
to assume a model for the halo profiles to connect the two mass def- 
initions. For this purpose, we used NFW profiles ( Navarro, Frenkj 
& White| 19971) with a scaling radius provided by the model of 
Bullock et al.| ( |200f| ) which relates to the virial mass of the halo 
through the concentration parameter c = ryir/fs. The M500 - Mytr 
conversion was obtained from the formulae provided in the Ap- 
pendix C of |Hu & Kravtsov| ( [2003| ), with the estimate of the over- 
density within the virial radius from |Kitayama & Suto| ( |1996| . 

The gas temperature within the dark matter halos is then com- 
puted using the M^qq-T relation of |Amaud, Pointecouteau & Pratt] 
(2005 ), i.e. a power low of slope a = 1.49, valid for clusters with 
T > 4 keV. To account for the steepening of the relation at lower 
temperatures suggested by their data, we used a higher slope of 
a - 1.9 below 3 keV. The lower halo mass limit is then fixed by 
imposing T > 0.5 keV. Bolometric luminosities are derived using 
the Lx-T relation of , Arnaud & Evrard] ( p'999 ). Self-similar evolu- 
tion is assumed for both these scaling laws. 

To account for the scatter observed in cluster properties, we 
encapsulate the dispersion of the M-T and L-T relations in the M-L 
relation, for the sake of simplicity. We assign the X-ray luminosity 
assuming a log-normal distribution. Following Stanek et al.] ( ]20061 », 
who measured crinM\L - 0.37, we use crinL\M ~ 0.37 X 1.59 ~ 0.6 
where 1.59 is the slope of their M - Lx relation. The total XMM- 
NEWTON EPIC count rate is estimated using the same spectral 
model as in section ]3.2] with a fixed hydrogen column density of 



2.6 X 10 cm~ (representative of our field), folded through the 
EPIC response matrices for the THIN filter in accordance with our 
observing mode. The selection function (Fig. [5]) is finally applied 
assuming, as in the previous section, a constant physical core radius 
of ISOkpc. 

In what follows, given the still modest size of the current clus- 
ter sample (some 30 objects over 5 deg^) we restrict most of the 
comparisons to qualitative visual ones. 



5.2 The redshift distribution 

Our model (using erg = 0.74, from WMAP3) predicts 6.2 Class 
1 clusters per square degree which, assuming Poisson noise alone 
(i.e.+l.l for the current 5.2 deg^), is compatible with our observed 
density of 5.4 deg~^ as given by the objects listed in Table [l] 

The observed redshift distribution of the CI sources, shown 
in Fig. [9] shows good overall agreement with the model expecta- 
tions, and suggests that we are crossing a void-like region within 
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Figure 9. Redshift distribution of the CI sample. The filled histogram 
shows the redshift distribution of our sample, while the solid line shows 
the expectations of our cosmological model (WMAP3: erg = 0.74 and self- 
similar evolution for the Lx -T relation). Expected fluctuations around the 
mean model density due to shot noise and (shot noise + cosmic variance) 
are respectively displayed with the solid and dotted error bars. The dashed 
line shows the expectations for a model with WMAP 1*^ year cosmological 
parameters (erg = 0.85) and a non-evolving Lx-T relation; combining this 
latter erg value with self-similar evolution would produce a normalisation 
even higher, by about 35%. 



Figure 10. The CI log(N)-log(S). The squares show the observed CI num- 
ber counts and the shaded region indicates the Icr error region. Overlaid 
are the predictions from our model: the dashed line for CI clusters, and the 
dotted line for the whole cluster population (i.e. all clusters above a given 
flux). Since only a fraction of cluster fluxes are detected in any survey, there 
is an important caveat in this comparison: our model delivers total clus- 
ter fluxes (i.e integrated to infinity) while our measurements are performed 
within 0.5 Mpc (which includes ~ 2/3 of the total flux for a cluster with 
Rc = 180 kpc and /3 = 2/3). Previously, Vikhlinin et al. (1998) claimed 
to measure total fluxes in an unbiased manner, while Rosati et al. (1998^ 
estimated that they had recovered at least 80-90% of the total flux. 



0.4 < z < 0.6. As is evident from the errors in Fig. [9] this ap- 
parent underdensity can be fully accounted for by statistical plus 
cosmic variance. We estimated the cosmic variance using the for- 
malism from jHu & Kravtsov| ( |2003| ), under the assumption that the 
surveyed volume is enclosed within a top-hat sphere; in the sensi- 
tivity regime of the XMM-LSS survey, the total uncertainty on the 
N{z) bins scales roughly as 2 Va^. 

Interestingly, the CI cluster density predicted assuming erg = 
0.85 (as inferred by |Spergel et al.|2003| from the V year WMAP 
data) is 14.6 deg~^ and appears to be much higher than required 
to match our observations. In this case, even using a non-evolving 
Lx - T, which is roughly our Icr lower bound from the previous 
section, our model predicts 10.7 CI clusters per deg^. Our data 
are 3.7cr (considering Poisson fluctuations) below this latter model. 
Even including cosmic variance, the difference is approximately 
2cr. This suggests that increasing the area by a factor of 2 would al- 
ready permit a stringent test of the value of erg indicated by WMAP. 

In this context, we examine in detail in Sec [63] some param- 
eters which play a role in the interpretation of the cluster number 
counts. 



5.3 The Flux and luminosity distribution 

As our selection function does not depend on flux alone, as as- 
sumed in previous generations of surveys, we have no direct esti- 
mate of the log(N)-log(S) relation. In this section, we provide the 
necessary information to enable the inter-comparison between our 
statistical distributions and those of past as well as future surveys. 

Our cosmological model can be both compared to the local 
log(N)-log(S), and used to predict an expected distribution for our 
CI sample . This is shown in Fig.[To| The raw log(N)-log(S) from 
our model is first compared with the estimates from |Rosati et al.| 



( |1998| ) and [Vikhlinin et al.||T998] ) based on flux limited samples, 
and said to be corrected for selection eflfects as well as for incom- 
pleteness. We find overall consistency, although our model suggests 
a somewhat higher total cluster density than the RDCS measure- 
ments above a flux of ~ 10"^"^ erg s~^cm~^. Furthermore, the pre- 
dicted flux distribution, once folded with the CI selection function, 
is consistent within the errors with the CI log(N)-log(S) relation 
derived from the present sample. This adds credibility to our full 
selection process and log(N)-log(S) modelling. 

Similarly, our sample is currently too small to derive the clus- 
ter luminosity function at any redshift, but this again can be pre- 
dicted by our model. We first compute as a cross-check the local 
luminosity function and compare it with the estimates from the RE- 
FLEX survey ( [Bohringer et al.|2002| , which constitutes the largest 
complete cluster sample currently available at low-redshift (Fig. 
[TTJ. The model shows, on average, good agreement with the mea- 
sured REFLEX luminosity function (which accounts for the survey 
sky coverage). It lies slightly below the REFLEX measurements 
for luminosities in the range 5 x 10'^^-5 x 10"^"^ erg s~\ behavior that 
is also present in the ROSAT Bright Cluster Sample (lEb eling et al 
[20001) as can be appreciated in the comparative plots of Mullis et al. 
( 2004 ). We further emphasize that, as is conspicuous in Fig.[TT] the 
inclusion in the model of scatter in the M -Lx relation is necessary 
in order to obtain predictions that are compatible with the REFLEX 
constraints at the high Lx end. Finally, we note that the prediction 
of our model at z = 0.8 (also displayed in Fig. [TTJ is very similar to 
the local one out to Lx = 10"^^ erg s~^ , and about a factor of 2 lower 
at Lx - 10"^^ erg s~^ (although still within the REFLEX error bars). 
Since the EMSS ( [Gioia et al[1990| ), the question of the evolution of 
the cluster luminosity function has been controversial, espec ially 
for luminosities above 5 x 10"^"^ erg s~^ [see a review by Mullis et[ 
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Table 4. Parameters of the mass-observable calibration from our cosmo- 
logical model, as defined by equations[To|and[TT] 



Parameter mass-flux relatior ^ mass - count-rate relatioi 



OBSo 

ao 
ai 
PO 
PI 



0.0329 
-37.68 
0.0645 
1.45 
0.248 



1.23 X 10^0 

-36.50 
0.0680 
1.41 
0.269 



^ For a total source flux in the [0.5-2] keV band in erg s ^cm ^ 
^ For a MOS1-FMOS2-FPN count-rate in the [0.5-2] keV band in cts 



[ar](f2004) including several ROSAT samples]. The most recent de- 
termination of the cluster luminosity function for Lx < 10^^ erg s~^ 
from the COSIVIOS-XIMIVI data shows no evolution ( [Finoguenov et| 
|al.|2007| ). Our calculations suggest that, interestingly, for a concor- 
dance cosmology, the observed luminosity function does not sig- 
nificantly evolve out to z ~ O.SforLx < 10"^^ ergs~\ while the 
cluster mass-function and scaling laws do. This can be regarded as 
a 'conspiracy' between cosmology and cluster physics. 




1042 10+^ 10+^ 10*5 

Luminosity [h-^ erg s"'' (0.1 - 2,4 keV)] 



Figure 11. The local luminosity function predicted by our model, compared 
with the REFLEX data points ( Boh ringer et al.|2002) . The dotted Hne is the 
best fit to the REFLEX data from the original article. The dashed line is 
the expectation from our model assuming no scatter in the M - Lx relation. 
The solid line shows our model (including scatter) yielding a much better 
agreement. The 3dot-dashed line is the prediction of our model for a redshift 
of 0.8. 



5.4 Mass-observable relations 

Last but not least, our model provides us with a tool for deriving 
heuristic mass-observable relations. For this purpose, we adopt the 
very general parametrisation introduced by Hu, ( 2Q03j : 

OBS = OBSoi — ] xe^^'\ (10) 

where Mj is the mass defined within some overdensity d, OBS 
is the observable of interest (flux, count-rate) and Mq and OBSq 
are reference values for both quantities. In practice, we used Mq = 
10^^ H'^Mq and Mj = M2oob i-^. the mass enclosed within the ra- 
dius delineating an overdensity of 200 with respect to the mean 
background density. IVlodel data points for a large range of masses, 
redshifts and fluxes (count-rates) were fitted with the above formula 
as illustrated on Fig. [12] We found that the functions p(z) and A(z) 
are sufficiently well described by functional of the form: 

/ p(z) = Po+PiZ 

\ A(z) = aoZ^^ ^ ^ 

The best fit values for the Mass vs Count-rate and Mass vs 
Flux relations are provided in Table |4] We do not give errors on 
the parameters as the overall accuracy of both fits with respect to 
the model data is better than 15% over the 0.05 < z < 1.2 and 
2 X 10^^ < M < 2 X 10^^ H'^Mq ranges. In practice, the intrin- 
sic dispersion of the M-L relation also needs taking into account. 
Assuming crinL\M ~ 0-6 (Sec. |5.1| l this translates into a dispersion 
of -45% -1-82% in the Flux- or Count-rate - Mass relation. Thus 
"perfect" flux measurements may yield, via this formula, mass ac- 
curacies of the order of -60% -hlOO%. Such an empirical relation 
has obvious useful practical usages. It should be noted however, 
that extrapolating the formulae in the current form above z > 1.2 
is not straightforward as a number of prominent low temperature 
lines {Oyu, Feyii, Feym) are redshifted below the 0.5 keV bound- 
ary, creating discontinuities that cannot be accounted for with the 
above simple functional. 





z = 0.10 


103 


z = 0.20 
z = 0.35 


102 


^^^^^ 2 = 0.65 - 


101 




10° 


^^^^^^^^^^^^^^^^^^^^^ 


10-1 




10-2 


1 1 , 11 ,~ 



10 100 1000 

M200b (1012 h-lMsun) 



Figure 12. Calibration of the [0.5-2] keV band total flux vs M^mh rela- 
tion. The exact correlations derived from our model are shown for several 
redshifts with plain lines. The overplotted dashed lines show the recovered 
values using the fit presented in Table|4] 

6 DISCUSSION 

We have shown that an accurate understanding of the sample se- 
lection effects is essential for a proper study of the evolution of the 
L-T relation of clusters of galaxies. Once these are allowed for, our 
data are compatible with self- similar evolution but cannot exclude 
no evolution at all. Subsequently, we have defined a cosmological 
cluster evolutionary model to which we applied the survey selec- 
tion function. It matches well our observations in terms of cluster 
number density, redshift and flux distribution, hence providing a 
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theoretical basis to further explore the ability of large XMM cluster 
surveys to constrain cosmology. 

In this respect, one of the critical issues is the role and the 
number of cluster parameters involved in the cluster scaling law 
evolution. As a first step, we have shown that for our best fitting 
model, it is possible to derive good mass-observable relations for 
the entire population of interest. We realistically explore below the 
impact of various parameters and hypotheses on the cluster scal- 
ing laws, on the selection function and, subsequently, on cluster 
number counts. In this way we hope to pave the way for future 
investigations. 



6.1 Cluster shape parameters 

Our selection function, as applied to our model, assumes a fixed 
core radius of = 180 kpc. Since Rc is related to the overall phys- 
ical size of the system, one naturally expects this value to depend 
on the cluster mass (or on any related parameter). This should ide- 
ally be included in the cosmological analysis, but we are currently 
lacking a well established Rc - T relation. The slope reported by 
|Ota & Mitsud a'('2004') is, for instance, much stronger than the one 
observed in the Sanders on & Ponman|(2003 ) sample. |Ota & Mit^ 
|suda| ( |2002| ) also revealed that, due to cool cores, the core radius 
distribution infered by performing a single yS-model analysis over 
the whole cluster population is actually double peaked and thus ill- 
defined. 

Furthermore, another known observational feature of X-ray 
groups and clusters is the tendency of the p parameter to drop with 
decreasing system temperature (see e.g. Jo nes & Forman|I999| or 
more recently [Osmond & Ponman 2004 ) but this correlation re- 
mains poorly defined, specially in the group regime. As our selec- 
tion function is mainly surface brightness limited, we expect the 
detection probability to depend on p as well as Rc and so, two ad- 
ditional scaling relations are needed. 

One could in principle derive them from the observed sample, 
but, given the generally small number of collected photons, it is not 
possible to resolve precisely the P - Rc degeneracy for most of our 
sources. Furthermore, such a study would be complicated by the 
fact that our selection process is dependent on these parameters, 
thus requiring a sophisticated self-consistent modelling similar to 
what was done in Section [43] for the L^-T relation. 

Given the small number of systems in the present sample, we 
postpone such a study to a future paper investigating the evolution 
of mean cluster/group profiles by stacking images (Alshino et al in 
prep.) 



6.2 The evolution of the Lx - T relation 

The constraints we were able to put on the evolution of the L-T 
relation (a = -0.07+j;;^J for the F{z) = E{z) X (1 + z)^ model) are 
not only limited by the relatively modest size of the sample and the 
uncertainties on the temperature measurements, but also depend on 
the assumed scatter of the relation (this parameter was fixed in the 
fit). 



6.2.1 Impact of the assumed scatter value 

In order to quantify the impact of the scatter in the determination 
of the evolution of the L-T relation, we ran a series of fits assign- 
ing different values to o-\^l\t- We find that the a index describing 



Table 5. Expected accuracy in the determination of the evolution of the 
Lx - T(z) oc E(z)il + zT relation for various sizes of XMM-LSS type sur- 
veys, estimated from simulations. The temperature measurements come 
either directly from the 10 ks XMM exposures, with corresponding accu- 
racy. Or, part of them are improved by subsequent deep XMM observa- 
tions, assuming an uncertainty of 10%; the last column gives the number 
of such clusters to undergo deep XMM pointings, era is the mean Icr error 
on a for a survey realisation. 



Area 


Temperature accuracy 


O-a 


Nre-obs 


5deg2 


from the survey 


0.59 




20 deg2 




0.28 




64 deg2 




0.15 




5deg2 


10% for the z > 0.8 clusters 


0.58 


2 


20 deg2 




0.25 


14 


64 deg2 




0.14 


41 


5deg2 


10% for the z > 0.5 clusters 


0.50 


9 


20 deg2 




0.22 


45 


64 deg2 




0.13 


139 



the deviation from self- similarity can be modelled by the following 
empirical formula: 

The uncertainty on crinriL of ~ 25% quoted by [Stanek et aL| 
( |2006| along with the Lx - T slope uncertainty (say between 2.5 
and 3) yield an uncertainty of ~ 20% for cti^lit (same reasoning 
as in Sec. |4.3.l] l. This translates to a conservative confidence range 
of [0.5 - 0.9] for o-\nL\T, which in turn gives fitted values of a in 
the range [0.6,-0.87]. The corresponding deviation is thus larger 
than the statistical errors from our fit, and precise constraints on 
cfin L\T are definitely required in order to place firm and rigorous 
constraints on F(z). Given that the scatter in the Lx - T relation 
could well be a function of redshift and mass, one will have to wait 
for large samples of high-z clusters with well monitored selection 
effects in order to undertake a fully self-consistent determination of 
the evolution of the slope and scatter of the Lx -T relation. 

6.2.2 Impact of the sample size and of the uncertainties on the 
temperature 

We have further investigated the extent to which the determination 
of the slope of the evolution of the Lx -T relation is conditioned 
by the size of the cluster sample and by the magnitude of the er- 
ror bars on the temperature (in comparison, errors on the luminos- 
ity are assumed to be negligible). We have thus simulated a large 
number of random cluster catalogues corresponding to 5, 20, and 
64 deg^ XMM-LSS type surveys. The realisations are drawn from 
our cosmological model, reproducing our observed «(z, T) distribu- 
tion along with selection effects; self- similarity is assumed for the 
Lx -T evolution. Each cluster temperature is assigned a mean er- 
ror estimated from the spectral fitting simulations presented in the 
appendix of Willis et al. ( 2005), for 10 ks XMM exposures. In ad- 
dition, we have considered the possibility of dedicated deep XMM 
follow-up observations providing a temperature accuracy of 10% 
for cluster subsamples of various sizes. For each sample, the slope 
was fitted following the method described in Sec |4.3.1| and |4.3.3| 
Results are summarized in Table |5] 

The simulations show that a significant improvement can be 
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reached by increasing the sample size as the accuracy on a scales 
roughly as the square root of the surveyed area. This is a very no- 
ticeable result, given that 2/3 of our clusters have no more than 
500 counts available for the spectral fit. The spectral accuracy is 
of the order of 20% below 2 keV, and 50% around 5 keV ( |Willis| 
|et al.|2005] ), the latter concerning mostly distant clusters. The sim- 
ulations further show that increasing the accuracy of the tempera- 
ture measurements does not yield a very significant improvement 
on the slope of the Lx - T relation - compared to the amount of 
time that would be necessary to obtain accurate temperatures for 
the 0.5 < z < 1 hot clusters. This is mostly due to the fact that the 
dispersion in the Lx-T relation itself is large. Quantitatively, apply- 
ing the yjarea rule of thumb on Fig.[7]3, shows that a 10 deg^ area 
(or any sample with a similar size and controlled selection eff'ects) 
has the potential to exclude the non evolution hypothesis. But even 
with 600 clusters (i.e. a 64 deg^ area, and 140 very well measured 
clusters above z > 0.5) it seems difficult to, for example, discrimi- 
nate between the two modified-entropy models of cluster evolution 
proposed by |Voit| ( [2005| . 



6.2.3 Working on samples close to the detection limit 

Fig. |8] shows that for the 0.2 < z < 0.4 range, most of our clus- 
ters are detected in a - T region where the selection effects 
are significant, as is to be expected, given that lower mass halos 
are much more numerous than massive ones. As a result, overlumi- 
nous clusters will be overrepresented in our sample. For the present 
data set, this is also true at any redshift, and this situation is easily 
understandable as a result of the combined eff'ects of (1) the high 
efficiency of Xamin close to the detection limit (it was designed 
to enable the construction of the largest possible controlled cluster 
sample), and (2) the steepness of the cluster mass function. Gener- 
ally speaking, given the contradictory former results obtained from 
data based on heterogeneous ROS AT- selected samples, our find- 
ings suggest that graphs such as Fig. |8] should be first carefully 
constructed, when studying the cluster Lx-T relation at any red- 
shift. 



6.3 Modelling the cluster number density 

While most of the ingredients linking cosmology to the expected 
cluster distribution are rather well established (mass functions, 
mean scaling relations for massive clusters) many details of the 
mass-observable issues are still pending. Among these, we would 
cite : the evolution of cluster scaling laws, the scaling laws in the 
group regime, and the role of the scatter in all relations. Our analy- 
sis has shown that our current data are compatible with self-similar 
evolution for the L-T relatiorQ In this section we thus focus on 
two other specific aspects, namely the slope of the M-T relation for 
the groups and the scatter in the M-L relation. We illustrate below 
the sensitivity of the cluster number counts to these quantities by 
means of our model. 
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Figure 13. Impact of the steepening of the M - T relation in the group 
regime on the modelling of the cluster number density. The full line corre- 
sponds to our fiducial model ( broken power law, Sec . |5.1) , giving 6.2 clus- 
ters per deg^. If we assume the |Vikhlinin et al| ( [2006) result (single slop e of 
1.5, dash line), the cluster density drops to 3.6 deg~^. The relation of |Ar-| 
|naud, Pointecouteau & Pratt ( 2005 1, motivated by a possible steepening of 
the mass function below 3.5 keV, (single slope of 1.7, dash-dot line) seems 
to be closer to our observations and yields a cluster density of 6.9 deg~^. 

6.3.1 The M-T relation in the group regime 

Our model assumes that the M - T relation for both massive clus- 
ters and groups can be modelled by a broken power law. Recent 
measurements | V ikhlinin et al| ( |2006] ) and |Arnaud, Pointecouteau &] 
Pratt ( 2005)] agree very well in the high mass regime, but yield con- 
tradictory results for groups. Vikhlinin et al (2006) claimed that the 
previously reported steepening of the relation (by e.g. |Finoguenov,| 
|Reiprich & Bohringer|2001 l) resulted from a bias toward low mass, 
due to an incorrect modelling of the density profile at large radius. 
However, their assertion relied on only two low-temperature sys- 
tems. On the other hand, Arnaud, Pointecouteau & Pratt| ( |2005| l 
found a steepening at low temperature based on a larger sample of 
groups, but their analysis assumes a NFW profile at large radius. 
As this point is currently still controversial, we adopted an inter- 
mediate solution: a steepening to 1.9 of the relation below 3 keV 
(note that a change in the Lx-T relation for groups would have 
similar impact). Fig, pj] illustrates how the various hypotheses on 
the shape of the M - T relation in the low-mass regime impact on 
the cluster number counts. Our assumed broken power-law model 
provides the best fit on the basis of visual inspection, but the dif- 
ference to the ' Arnaud, Pointecouteau & Pratt] ( |2005i ) model is not 
statistically significant with the current sample size. 

6.3.2 Impact of the scatter 

We discussed above how the scatter in the Lx-T relation influ- 
ences the determination of the evolution of the relation from the 



^ In this respect we stress that combining a fully self-similar M-T relation 
with 2i Lx - T relation whose slope is not the one predicted by the self- 
similar model (only gravitational physics) implies that the Lx-T and M - 
Lx relations cannot simultaneously evolve in a self-similar way. The point 
seems to have been overlooked in the cluster detailed studies so far 
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Figure 14. Impact of scatter in the M - Lx relation on the predicted CI 
redshift distribution. For given scatter values, the predicted total number of 
clusters per deg^ are respectively : 3.6 (no scatter), 5.1 (0.3), 8.4 (0.9); our 
model gives 6.2 (0.6) . 



observations (Sec. |6.2.1| ) . In our model, all sources of dispersion 
are encapsulated in the M-Lx relation, since the determination and 
understanding of dispersion in Lx-T and M-T, and thus M-Lx, 
and of their possible evolution, is still in its infancy. In this respect, 
we note that measurements of scatter other than those of St ane k et 
[aT (|2006|) are availab le: lower values for cti^lim have been found 
by Zhang et al.|2006| and |Zhang et al. 12007] (respectively 0.41 and 
0.33), but they are dominated by higher temperature systems which 
are known to have more regular properties. In terms of the accuracy 
with which the scatter is determined, the typical errors estimated by 
[Stanek et al.| ( |2006| ) lead to an uncertainty in o-\nL\M of the order of 
20%. 

In Fig. [14] we show the impact on the CI redshift distribu- 
tion of neglecting the scatter in M - Lx, as well as the influence 
of plausible values of scatter. If the luminosities are symmetrically 
scattered around the expected value for a given mass, the net effect 
on our cluster counts (approximately the number of clusters de- 
tected above a given luminosity at a given redshift) is an increase of 
the detected cluster density since, due to the steepness of the mass 
function, substantially more low-mass clusters become detectable 
than high-mass cluster undetectable. 



6.4 Cosmology and duster survey self-calibration 

We have demonstrated that taking into account survey selection 
effects is not only critical for the proper modelling of the cluster 
number counts (hence cosmology) but also for the determination 
of the cluster evolutionary scaling laws, which in turn, also impact 
on dn/dz (Fig.|9]). Further, we explored the impact on dn/dz of the 
slope of the M - r relation (Fig. [13]) and of the dispersion in the 
scaling laws (Fig. [14]). This illuminated, by means of real data, the 
magnitude of some of the degeneracies between cosmology and 
cluster evolution. However, the number of cosmological free pa- 
rameters that we have explored in the present study is very lim- 
ited (basically, only erg, Fig.|9]) and a more general analysis should, 
among others, include constraints on the Dark Energy parameter w. 
For instance, in Fig. [15] we exhibit the predictions of our model for 



different values of w and a (evolution parameter of Lx-T). A clear 
degeneracy is apparent, adding to those already reviewed. 

In order to cope with this critical issue, it has been suggested 
that X-ray (or SZ) cluster number counts, in principle, have the 
ability to self-calibrate even an evolving mass-observable relation, 
because the number density of clusters as a function of redshift has 
a fixed functional form given by c osmological simulations (|Levine,| 
[Schulz & White||2002|[H"u||2003|[Majirmdar & Mohr|2003j ). This 
suggests the possibility of simultaneously solving for the cosmol- 
ogy and fitting an ad hoc parametrized mass-observable relation, 
provided that a sample of several thousands of clusters is available. 

While our - still somewhat phenomenological - best fitting 
model, allowed us to derive flux- and count-rate - mass relations, 
several caveats still have to be pointed out. We have noticed that the 
flux-mass relation is very sensitive to the details of the X-ray emit- 
ting plasma (e.g. emission lines) and also, obviously, to the way 
fluxes are estimated in shallow XMM exposures; in earlier studies, 
even the temperature-dependence of the bremsstrahlung continuum 
was ignored (e.g. |Majumdar & Mohr 2003 ). Even neglecting mea- 
surement uncertainties, our mass - observable mapping exhibit Icr 
uncertainties in the range [-60%,+ 100%] coming solely from the 
dispersion of the M - Lx relation. 

In this respect, while vastly increasing the number of clus- 
ters may be useful, additional independent observations on cluster 
evolutionary physics will constitute a crucial input. This will im- 
prove the modelling of the mass-observable relations, and thus the 
constraints on the mass function. In particular, deeper XMM point- 
ings will provide accurate temperatures and, consequently, a better 
estimate of the dispersion in the L-T relation. Weak-lensing, using 
the CFHTLS images, will provide independent mass measurements 
(Berge et al in prep.). 

In the future, the XMM-LSS field will be covered by Sunyaev- 
Zel'dovich (SZ) observations [OCRA (One Centimeter Receiver 
Array), AMiBA (Array for Microwave Background Anisotropy) 
and APEX (Atacama Pathfinder Experiment)]. Since the SZ decre- 
ment is an integral of the cluster pressure (rieTe) along the line of 
sight, and is independent of z, it provides an especially interesting 
complement to the X-ray emissivity (scaling as nlT^^^ for moder- 
ately hot clusters). Comparing our limiting mass as a function of 
redshift with those of these SZ surveys is all the more informa- 
tive as we introduced a new class-oriented selection to define the 
XMM-LSS cluster sample. This is done, as a test case, in Fig. [T6]fo r 
the XMM-LSS and the APEX-SZ surveys [Schwan et al. (2003]), 
[Dobbs et al.| ( |20061 )]. The SZ limit has been obtained under the as- 
sumption that APEX will observe at two frequencies (150 and 220 
GHz) with a V beam and 10 yt/K/beam white noise in each chan- 
nel. We have included the effect of contamination by primary Cos- 
mic Microwave Background (CMB) anisotropics. Given the spatial 
spectrum of the CMB, the instrumental noise and cluster profiles, 
we compute the expected SZ flux limit as a function of cluster size 
using matched filters. We then convert from SZ fluxes and sizes to 
masses and redshifts, using the same cosmology and scaling laws 
as in our X-ray modelling. Details of the computation are given in 
IMelin et al.|pbo5] ). 

The XMM-LSS mass limit of 10^"^ around z = 1 for the 
CI selection, appears comparable to the prediction for the deepest 
SZ surveys to date. In addition, the XMM-LSS uncovers more low- 
mass systems below z < 0.5. According to our current model, the 
APEX-SZ cluster density detected at 3cr is 3.8 deg~^ out to z = 1, 
to be compared with 5.4 deg~^ for the observed CI X-ray selection. 
In a forthcoming paper (Melin et al in prep.) we explore the cos- 
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Figure 15. Dark energy vs scaling laws for our modelled cluster reshift 
distribution. The solid line shows the prediction for the jSpergel et al.|{2006) 
cosmology and self-similar evolution of the scaling laws. Assuming instead 
F(z) = E(z)x(l+z)'^-^^ for the Lx-T evolution, which is our Icr higher bound 
from section [43] yields the dotted-dashed line. As a comparison, the dashed 
line shows the distribution obtained by keeping self-similar evolution but 
changing the dark energy equation of state from -1 to -0.6. 
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Figure 16. Comparison of the XMM-LSS and APEX-SZ sensitivity in 
terms of limiting mass. The red lines show various probability detection 
thresholds for the CI clusters. The blue lines are the predictions for the 
lOjdK APEX survey (see text). The X-ray and SZ curves make use of the 
best fitting model to our current data set (Sec. |5.1) . 

mological constraints expected from a few hundred clusters, whose 
masses are determined by a joint X-ray/SZ/weak-lensing analysis. 



7 SUMMARY AND CONCLUSION 

• This work presents a well-controlled XMM cluster sample 
over 5 deg^ with a density of 5.4 deg~^ at medium X-ray sensitivity 
(down to ~ 10"^"^ ergs~^cm~^ for the extended sources in ques- 
tion). A complementary sample of about the same size, but with 
less well defined selection criteria is in preparation (Adami et al). 



We provide positions, redshifts, fluxes, luminosities, temperatures 
and masses along with X-ray and optical images. 

• The selection is based on well-defined criteria pertaining to 
spatial extent properties, and are similar to a surface brightness 
(rather than flux) limit. 

• The resulting cluster redshift distribution extends out to z = 
1.2 and peaks at z ~ 0.3. Half of the clusters have temperatures in 
the range 1-2 keV, occupying the 0.2 < z < 0.4 interval. This inter- 
mediate population, the building blocks of the present day clusters, 
is systematically unveiled by the XMM-LSS survey. It is also the 
first time that a wide- area blind X-ray survey has provided reliable 
cluster temperature measurements. 

• We demonstrate that taking into account the sample selection 
eff'ects is vital for a proper determination of the evolution of the 
Lx-T relation. This is due in large part to the steepness of the 
cluster mass function, which results in sources clustering close to 
the detection limit, and of these the overluminous systems are pref- 
erentially detected. This may explain the often contradictory re- 
sults obtained with heterogeneous samples selected form ROSAT 
surveys. 

• Modelling the selection eff'ects carefully, our data appear to 
point to approximately self-similar evolution. A sample twice as 
large should allow exclusion of the no-evolution hypothesis, at bet- 
ter than the Icr level. Using extensive simulations, we find that in- 
creasing the sample size is more efficient than increasing the ac- 
curacy of the temperature measurements for constraining the evo- 
lution; this is due to the large intrinsic dispersion of the Lx-T 
relation. 

• Our results suggest a higher normalisation of the cluster 
log(N)-log(S) at faint fluxes than previously obtained by deep 
ROSAT surveys. 

• Cluster masses were estimated from surface brightness profile 
fitting, assuming hydrostatic equilibrium. Our sample follows sur- 
prisingly well (given the very modest XMM exposures of ~ 10 ks 
and our rudimentary modelling of the mass profile) the local Lx-M 
relation, when evolved self- similarly. 

• Self-consistent cosmological modelling of the cluster popu- 
lation, convolved with the accurately determined survey selection 
function, confirms that the properties of our current data set are 
compatible with the concordance cosmology along with cluster 
self- similar evolution. 

• This model then allowed us to investigate several degeneracies 
arising from cluster physics, regarding the cosmological interpreta- 
tion of the number counts. In particular, we stress the influence of 
the, still poorly determined, scatter in the cluster scaling laws. 

• The present study led us to investigate many of the issues 
raised when attempting to perform precision cosmology using real 
cluster data, including the need for precise selection eff'ects and a 
realistic error budget. Thanks to the temperature information ob- 
tained from the detected photons for all CI clusters, this work con- 
stitutes one of the first attempts to break, in a self-consistent way, 
the degeneracy between cosmology and cluster evolution, in the 
analysis of the cluster number counts. The next step will be to ap- 
ply cosmological modelling to the full XMM-LSS area (10 deg^ 
- to be completed in 2007) adding input from the combined weak 
lensing and SZ survey as well as from the 3 -dimensional cluster 
distribution. 

• All data presented in this paper - cluster images taken 
at X-ray and optical wavebands in addition to detailed re- 
sults for the spectral and spatial analyses - are available 
in electronic form at the XMM-LSS cluster online database: 
http : //13sdb . in2p3 . fr : 8080/1 3 sdb/ 
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APPENDIX A: NOTES ON SPATIAL FITTING AND 
PHOTOMETRIC ACCURACY 

In this appendix, we test the accuracy of our luminosity measure- 
ment procedure, running the spatial fitting algorithm on simulated 
images. 

For this purpose, we created 10"^ 5 XMM pointing simulations 
with vignetted and particle background levels set to the mean ob- 
served values of Read & Ponman, ( ,2003^ , and point sources dis- 



tributed according to the log(N)-log(S) of |Moretti etaL] ( |2003| ). We 
then randomly included in these simulations extended sources with 
p -213 and several Rc and count rate values (respectively from 10 
to 50''and 0.005 to 0.1 cts s"0. A /5-profile was fitted to all the sim- 
ulated sources that were detected by the first pass of our pipeline 
(see jPacaud et al.|2006| for details) and for which the spatial fitting 
procedure detects at least three radial bins above 3cr. 

The simulations showed that that below 400 detected counts 
in the spatial fitting radius, the degeneracy between p and Rc is 
large, but that the integrated flux (within the fitting radius) is well 
modelled by the best fitting (fi,Rc) combination. Results of the sim- 
ulations are presented in Fig. |Al| 

The measurements are almost unbiased: less than a few per- 
cent in general, and up to ~ 10% for very faint sources (100- 
200 counts). This small off'set can be interpreted as a weak Edding- 
ton bias, as it is of the order of Poisson fluctuations, and proba- 
bly results from our requirement of recovering at least 3 significant 
bins (which means retaining only the clusters that appear brighter 
among the faintest ones). This shouldn't be an issue for the CI 
clusters as they are generally brighter (see Table [T]). The scatter 
around the mean value is very low for bright sources and mildly 
increases (up to ~ 15%) for the faintest ones. Estimating the lumi- 
nosity within a diff'erent radius than the fitting one does not result 
in significant diff'erences (although the scatter increases) as long as 
the profile is not extrapolated too far out, which is the case for our 
sample. Finally, the mean estimated errors are of the same order as 
the dispersion of the best fit value, assessing our confidence inter- 
vals. 

Other sources of error have already been considered by |Pierre| 
|et al.| ( [2006] ) appendix A: the impact of neglecting the errors on the 
temperature (which aff'ects our estimate of R500) was shown to be 
smaller than a few percent, while the contribution from undetected 
weak AGNs is probably lower than the percent level (which is com- 
forted by the present simulations in which AGNs are included). 

We thus conclude that the possible systematics generated with 
our fitting procedure are rather low, especially compared to the sta- 
tistical errors. 

APPENDIX B: INDIVIDUAL CHARACTERISTICS OF 
THE CI GALAXY CLUSTERS 

Notes on individual sources: 

• XLSSC-028: 

The velocity information obtained to date on this object does not 
allow us to firmly conclude on its redshift. The two bright galaxies 
on which the X-ray emission is centred have a redshift of 0.08; their 
spectra are typical of elliptical galaxies, without emission lines. No 
other galaxies with this redshift have been measured in the field, 
but a number of z=0.3 objects are found within 500 kpc of the X- 
ray centroid. The measured X-ray temperature is 1.5 keV for z=0.3 
and, the luminosity is found to be \.6e43 ergs~\ which puts the 
object close to the observed L^-T relation. At z = 0.08, the source 
is about a factor of 17 fainter, and corresponds to a temperature of 
0.75 keV which put it also exactly on the L^-T relation. 

• XLSSC-018: 

The first measurement of temperature of this system (~ 2.7^q| keV, 
published in |Willis et al.| ( |2005] )) was one of the hottest in our sam- 
ple at a redshift around 0.3. Yet, the cluster doesn't seem mas- 
sive, neither in the X-ray nor in the optical. The emission in the 
[2-10] keV appears to be significantly displaced from the low en- 
ergy centroid (Fig. |Bl| ), and is possibly associated with an optical 
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Figure Bl. Gaussian smoothed X-ray emission from source XLSSC-018 at z=0.32 in bands [0.5-2] keV (left) and [2-10] keV (right). The emission in the soft 
band seems slightly offset from the cluster cD galaxy. In the hard band, the offset becomes much larger. 



counterpart. This suggests that a fraction of the emission could be 
due to an AGN. Excluding the probable contaminated region from 
the spectral fit lowers by a factor of two the number of available 
photons, but surprisingly leads to much tighter error bars around a 
temperature of 2.0 keV (as reported in Tabl^. We regard this as 
strong evidence that the cluster's X-ray emission suffers AGN con- 
tamination, and consequently discard this source from the Lx - T 
analysis. 

• XLSSC-006: 

This source was already studied in [Willis et al.| ( |2005| ), and 
assigned a luminosity of (4.5 + 0.3) x 10^^ ergs~^ which is 
incompatible with our estimated (6.0 + 0.2) x 10^^ ergs"^ In that 
paper, the luminosity was derived from the Xspec spectral fitting 
measurement, to which a correction factor was applied in order to 
extrapolate to rsoo- We interpret the discrepancy as resulting from 
the combination of two effects. Firstly, Xspec uses a background 
estimate comming from an annulus surrounding the source, while 
in our spatial fitting procedures, we model it over the whole 
pointing. As a consequence, the mean background level used 
within Xspec actually happened to be higher than our fitted value 
due to residual contamination of the bright XLSSC-006 cluster 
in this local annulus. Additionaly, the observed source count-rate 
in the [0.5-2] keV band (on which our spatial fitting luminosity 
estimate is based), is found to be higher within Xspec than the one 
infered from the best spectral fitting model. 

Information on other CI sample clusters is already published 
by|Valtchanov et al.|p004| ), , Willis et al.,(2005j and Pierre et al.| 
( [20061 l 



APPENDIX C: THE CI NEARBY GALAXIES 

These sources were discarded from the CI galaxy cluster sample 
based on an obvious link between the main X-ray emission, and the 



Table CI. The nearby galaxy sample. XMM pointing identifiers refer to the 
XMM-LSS internal naming convention as described in Pierre e t al.| ( [sub^ 
[mitted) ; the location of each pointing on the sky is shown in figure [T] 



Source name Pointing R.A. Dec. 

(J2000) (J2000) 



XLSS J022528.7-040041 


B03 


36.3699 


-4.0115 


XLSS J022251.4-031151 


Bll 


35.7146 


-3.1975 


XLSS J022659.2-043529 


G06 


36.7469 


-4.5916 


XLSS J022430.4-043617 


G08 


36.1268 


-4.6048 


XLSS J022617.6-050443 


G16 


36.5735 


-5.0788 



presence of a nearby galaxy on the same line-of- sight. A significant 
fraction of the total X-ray emission can however originate from 
another object in the field. The source list is given in Table |C1| and 
details regarding the origin of the X-ray emission are given below: 

• XLSS J022528.7-040041: 

The nearby galaxy lying a the centre of the X-ray isophotes is 
alreday known as AP]V[UKS(BJ) B022258.83-041412.5. The X- 
ray emission from this galaxy certainly originates from its inter- 
action with two satellite companions, as conspicious in the images 
ofFig.jcT] 

• XLSS J022251.4-031151: 

This source is likely to result from the confusion of a point source 
with the nearby elliptical galaxy APMUKS(BJ) B022018.98- 
032531.1. It is associated with a radio source detected i mTasse eti 
|aL H2006l ). 

• XLSS J022659.2-043529: 

The nearby galaxy lying at the centre of the X-ray isophotes is al- 
ready known as IMCG -01-07-011. We cannot exclude it to be the 
dominant galaxy of a poor nearby group, as its aspect is reminiscent 
of the compact group XSLC-01 1 at z 0.05. 
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• XLSS J022430.4-043617: 

This source is likely to result from the confusion of a point source 
with a nearby spiral galaxy. The galaxy is already known as 
6dF J0224300-043614, and has a redshift of z=0.06916. 

• XLSS J022617.6-050443: 

The nearby galaxy lying a the centre of the X-ray isophotes (AP- 
MUKS( BJ) B022347 .26-051811.5) is detected as a weak radio 
source in |Tasse et al.1 ( |2006 ). As for, XLSS J022659.2-043529, we 
cannot exclude the possibilty that it is the dominant galaxy of a 
small group. Another option would be that it is part of the recently 
identified XBONG class (see e.g. |Yuan & Narayan|2004l ). 
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Figure B2. Images of the CI clusters (sorted by right ascension). Left: X-ray/I band overlay of the central 1' . Right: true colour image with X ray contours 
and measured redshift in the central 1.5 Mpc. (a): XLSSC-039. (b): XLSSC-023. 
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ARC i^iNuies 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-006. (b): XLSSC-040. 
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AftC MINUTES 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-036. (b): XLSSC-047. 
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ARC w\m\E5 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-048. (b): XLSSC-035. 
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ARC I^INUICS 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-028. (b): XLSSC-049. 



© 2007 RAS, MNRAS 000,[l](38] 



The XMM-LSS Class 1 sample over the initial 5 deg^ 




Figure B2 - continued Images of the CI clusters, (a) XLSSC-018. (b): XLSSC-029. 
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-J -2 -t 1 ? 5 

ARC t/inuTE£ 




ARC 4^INU1C£ 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-044. (b): XLSSC-021. 
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ARC 4^INU1C£ 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-001. (b): XLSSC-008. 
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ARC ^INUHK 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-025. (b): XLSSC-041. 
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ARC t^lNUlCS 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-002. (b): XLSSC-050. 
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Figure B2-co^?mMg J Images of theCl clusters, (a) XLSSC-051. (b): XLSSC-011. 
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ARC wmtui 

Figure B2 - continued Images of the CI clusters, (a) XLSSC-052. (b): XLSSC-005. 
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Figure B2 - continued Images of the CI clusters, (a) XLSSC-010. (b): XLSSC-013. 
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Figure B2 - continued Images of the CI clusters, (a) XLSSC-003. (b): XLSSC-022. 
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Figure B2 - continued Images of the CI clusters: XLSSC-027. 
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Figure CI - continued Images of the CI nearby galaxies, (a) XLSS J02225 1.4-03 1151. (b): XLSS J022659.2-043529. 
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Figure CI - continued Images of the CI nearby galaxies, (a) XLSS J022430.4-043617. (b): XLSS J022617.6-050443. 
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